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gions where the expected upward-going muon fluxes are within the expected reach 
of AMANDA II. 
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1. Introduction 

There is now compelling evidence for a non-baryonic cold dark matter (DM) compo- 
nent in the universe [Q]. In the Minimal Supersymmetric Standard Model (MSSM) 
with R parity, the lightest supersymmetric particle (LSP) is stable and becomes a 
good candidate for cold dark matter in the universe 0. The LSP is often the lightest 
neutralino which is the admixture of Bino, Wino and Higgsinos in the MSSM. In this 
case, the neutralino DM in our galactic halo might be detected via its elastic scatter- 
ing with terrestrial nuclear targets @, [|. In fact, the DAMA Collaboration || even 
claimed an evidence for DM. However, the CDMS II experiment has reported the 
upper limit on the DM scattering cross section, which is not compatible with the 
results of the DAMA experiment. There are several experiments going on searching 
for the DM scattering at the level of of™ ~ 1CT 7 pb or less. In the most widely 
studied minimal supergravity (mSUGRA) scenario (or the constrained MSSM), the 
spin-independent neutralino-proton scattering cross section a xp turns out very small 
(< 10~ 8 pb). However there is no solid theoretical rationale for the minimal super- 
gravity scenario, and it is important to calculate the possible maximal values for 
<j xp in general supergravity scenarios beyond the mSUGRA scenario. And it is very 
important to impose all the relevant constraints from various experiments in order 
not to overestimate the cross section. Some important constraints include the lower 
bounds on the Higgs and SUSY particle masses, B —>■ Xg^ branching ratio, the muon 
(g — 2) M , etc.. One may also take some theoretical consideration on the absence of 
the color-chrage breaking minima or the directions unbounded from below, etc.. 
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In a previous work , we pointed out that there is a strong correlation between 
the spin independent neutralino-proton scattering cross section a xp and the branching 
ratio for B s — ► decay || |9], |T0], [TT|, [L2[] within mSUGRA and its extensions. 

The origin of this correlation resides in the dependence of both observables on tan /3 
and the neutral Higgs boson masses; both observables increase for large tan/3 and 
low Higgs masses. In particular, we have shown that the current upper limit on 
B(B S — ■> excludes substantial parameter space where the DM scattering cross 

section is within the CDMS sensitivity region (see also |~3| for a detailed analysis). 



In this work, we extend our previous study to the indirect detection of neutralino 
DM with neutrino telescope through upward-going muon flux, and its correlation 
with B(B S — > fi + fi~). The energetic neutrino(-induced muon) flux from neutralino 
DM annihilation in the sun and the earth is one of the promising signals in the 



indirect detection of neutralino DM [13]. Neutralino DM particles in the halo can 
be captured by the sun or by the earth, when their velocities drop below escape 
velocities via their elastic scattering with matter in the sun or earth. Then they 
will accumulate in the core of the sun and the earth and will eventually annihilate 
into ordinary SM particles. Among the annihilaion products, neutrinos can pass 
through the sun and the earth, and then could be detected in neutrino telescopes 
through their conversion to muons via charged- current scattering with neuclei near 



the detectors. Baksan pf, MACRO fl§, Super-K |TJ and AMANDA [|T|] released 



upper limits on the upward-going muon flux. There are also planned or proposed 



neutrino telescopes such as ANTARES [19], IceCube [20] and NESTOR [EJ etc.. 



An important point of the indirect detection of neutralino DM with neutrino 
telescopes is that the neutrino flux strongly depends on the capture rate of neutralino 
by the sun or the earth, which in turn depends on neutralino-nucleon scattering 
cross sections. Therefore we expect some correlation between the neutrino flux and 
B(B S — > /i + /i _ ), which is similar to the strong correlation between a xp and B(B S — > 
/i + yU~) as discussed in Ref. @. Indeed, we will show that the current upper limit 
of B(B S -> fi + fi-) < 4.1 x 1(T 7 (90 % CL) [pj, || puts strong constraints on 
the upward-going muon flux in the supersymmetric models which give rather large 
spin-independent neutralino-proton scattering cross section. 

This paper is organized as following. In Sec. 2, we give a brief review on 
the indirect detection of the DM through the upward-going muon flux. In Sec. 3, we 
consider the upward-going muon fluxes in some supergravity scenarios and illuminate 
our point that B s — ► fi + n~ branching ratio plays an important role. In Sec. 4, we 
summarize the results. 

2. Indirect detection through the upward-going muon flux 

As we mentioned in the introdecution, the observation of energetic neutrinos from the 
sun and/or the earth would provide convincing evidence of the existence of neutralino 
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dark matter in galactic halo [|J . The flux of energetic neutrinos from neutralino anni- 
hilation in the sun or the earth is proportional to the rate of neutralino annihilation 
in the sun or in the earth and the energy spectrum of neutrinos from the annihilation. 
The time evolution of the number of neutralino, N in the sun (or in the earth) is 
given by 

N = C-C A N 2 (2.1) 

where C is the capture rate of neutralino by the sun or the earth and C A is the total 
annihilation cross section times relative velocity per volume. From Eq.(2.1), we find 
that the present annihilation rate is 

T A = \c A N 2 = 1 -C%mh 2 {^/CC~ A t ) (2.2) 



where to — 4.5 Gyr is the age of the solar system. For ^JCCa to <SC 1, the annihilation 
rate is Ya ~ \C 2 Ca^ and less than its maximal value. But, for \JCCa to ^ 1; the 
neutralino density reach equilibrium and the annihilation rate is Fa ~ \C. Therefore, 
when accretion is efficient, the annihilation rate depends on the capture rate C, but 
not on the annihilation cross section. 

In turn, the capture rate C strongly depends on the elastic scattering cross 
section of neutralino with matter in the sun and the earth. The capture rate for 
the earth primarily depends on the spin-independent DM scattering cross section. 
For the capture rate in the sun, however, both spin-independent and spin-dependent 
scattering cross section can be important and the significance of each contribution 
depends on the specific SUSY scenarios. 

In MSSM, t-channel Higgs boson and s-channel squark exchange processes con- 
tribute to the spin-independent (scalar) scattering between neutralino and quarks. 
In many case, dominant contribution to the scalar cross section comes from the 
Higgs exchange process, which increases for large tan j3 and small Higgs masses and 
also if neutralino is a mixed gaugino-Higgsino state. On the other hand, for the 
spin-dependent cross section, t-chennel Z boson and s-chennel squark exchange pro- 
cesses contribute. Usually Z exchange contribution dominates, which is sensitive 
to Higgsino components of LSP, but largely independent of tan (3. Note that if the 
Higgsino component of the LSP increases, then both the spin-independent and the 
spin-dependent scattering cross sections will be enhanced, as shown below in the 
nonuniversal Higgs mass parameter case. 

The capture rate C also depends on the local density of neutralino, p x and 
the neutralino velocity dispersion in the halo, v etc. For our numerical calculation, 
we use the code DARKSUSY [^4] and fix v = 270 km/s. For the local density of 
neutralino we fix p x = 0.3 GeV/cm 3 if fl x h 2 > 0.025, while performing a rescaling 
of the density as p x — > p x (Q x h 2 /0.025) if £l x h 2 < 0.025. Here fl x is the neutralino 
relic density in units of the critical density and h is the present Hubble constant in 
units of 100 km s -1 Mpc -1 . 
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3. Upward-going muon flux in SUSY models 



3.1 mSUGRA 

In mSUGRA model, we assume a universal SUSY breaking scalar mass m, a universal 
gaugino mass M and a universal trilinear coupling A at GUT scale wigut — 2 x 10 16 
GeV. We also require that electroweak symmetry break radiatively and then the 
Higgsino mass parameter /i is determined by the condition: 

„ m 2 TT — m 2 TT tan 2 3 1 „ 

where tan/5 is the ratio of two Higgs vacuum expectation values and m 2 Hu1 m 2 Hd are 
the soft breaking Higgs masses-squared. With the above mSUGRA assumptions, 
\n\ is usually large so that the lightest neutralino is bino-like and the pseudo-scalar 
Higgs mass uia is rather large. 




(a) $ M (sun) (b) $ M (earth) 

Figure 1: The muon flux from the sun and the earth vs. m x in mSUGRA scenarios. The 
three branches correspond to tan/3 =10,35 and 50 respectively (from bottom to top). The 
current upper limit from SUPER-K and the expected reach of AMANDA II are also 
illustrated. 

In Fig.s [l] (a) and (b), we show the allowed ranges of the upward-going muon 
fluxes from the sun and the earth respectively as functions of the LSP mass. The 
three branches correspond to tan/5 = 10,35 and 50 cases (from the bottom to the 
top), respectively. Here, we took A = and \i > (motivated by the muon (g — 2) M 
experiment) and varied m and M up to 1 TeV. We have imposed the experimental 
bounds for the Higgs and sparticle masses and for b — > branching ratio. We 
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also required that the lightest neutralino is LSP. For opposite sign of /x, the muon 
flux could be smaller. But we are interested in the possible maximal values, and we 
consider the positive /i case only in this work. 

In our scan, the muon flux from the sun reaches up to ~ 20 km~ 2 yr~ l if the neu- 
tralino LSP is light enough (m x ~ 100 GeV). In the small m x region, the neutralino 
density in the sun can reach (near) equilibrium so that the muon flux is more or less 
determined by the capture rate of neutralino by the sun. And in turn, the capture 
rate in the sun is determined primarily by the spin-dependent scattering cross sec- 
tion, (though the contribution from the spin-independent scattering cross section to 
the capture rate can be comparable to the one from the spin-dependent cross section 
for very large tan/3 cases) As we already mentioned in Sec. 2, the spin-dependent 
scattering cross section is largely independent of tan (3. Therefore the upward-going 
muon flux from the sun in the small m x region gives similar values for the three 
choices of tan f3 values, as one can check from the figure. 

In large m x region, the neutralino density is usually far from equilibrium since 
the elastic scattering cross section of the DM with ordinary matter in the sun be- 
comes smaller, and the neutralino annhilation cross section becomes important for 
the prediction of the muon flux. An important process in this case is the neutralino 
pair annihilation into bb through s-channel pseudo-scalar Higgs exchange diagram, 
which is strongly enhanced for large tan/? ||26|| . From the Fig. |l] (a), one can notice 
a clear dependence of the muon flux from the sun on tan j3 in the large m x region. 

For the muon flux from the earth, the neutralino density is far less than the 
equilibrium values and both the capture rate and the annhilation rate are important 
for the calculation of the muon flux. The resulting flux is much below the one from the 
sun. The maximal value of the muon flux from the earth is about 3 x 10~ 5 km~ 2 yr~ l , 
which is far below the SUPER-K and AMANDA II sensitivity regions. 

Note that there is no further constraint from the B s — > bound for the 

mSUGRA case, once we impose the constraints from the lower bounds for Higgs 
boson and SUSY particle masses and the B — ► A s 7 branching ratio, as discussed in 
Ref. 0. 

3.2 Non-universal Higgs model (NUHM) 

In the previous subsection, we have shown that the mSUGRA assumption predicts 
the muon fluxes from the sun and the earth that are far below the sensitivity region 
of the current experiments. This is mainly because the lightest neutralino is bino- 
like and the pseudo-scalar Higgs mass is large in the scanned region of mSUGRA 
scenario. Larger muon flux from the sun and the earth can be obtained if we relax 
the universal boundary condition at GUT scale. 

In this subsection we consider the non-universal Higgs model, in which the as- 
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sumption of universal soft scalar masses are relaxed for soft Higgs masses, as follows: 



ni 



2 Hu =m 2 (l + 5 Hu ) 



ni 



m 2 (1 + 8 Hh 



(3.2) 



whereas other scalar masses still have a universal mass m at GUT scale. Here <5's 
are parameters with < 0(1). 
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Figure 2: a^ ar vs. a s x p p m /(2m x /GeV) in NUHM (black) and in mSUGRA (green) for 
(a) tan /? = 35 and (b) tan (3 = 50. 

As an optimal choice for enhancing the muon flux from the sun and the earth, 
we take the numerical values of 5's as 5u d = — 1 and 5h u = 1- For the postive 5h u , H 
becomes lower and the Higgsino component in the neutralino LSP increases so that 



a xp is enhanced, as discussed in Ref. |27], |28|]. The change of also has an impact 
on the Higgs masses because 



ml + m 2 + 2fi 2 ~ m 2 + /i 2 - M 2 z /2 



at weak scale. For the negative SH d , fn>A and tuh become further lower. As the result, 
both spin-independent scattering cross sections <j s x c £ lar and spin-dependent one cr^ m 
are enhanced compared to mSUGRA case. 

In Fig. |, we present a s x c p alar vs. a%£ n /{2m x /GeV) in the NUHM (black points) 
and mSUGRA scenario (green points) for (a) tan /3 = 35 and (b) tan /3 = 50 respec- 
tively. The dashed straight line in the figure indicates the region in which the two 
contributions to the capture rate are similar to each other. We observe that both o s x v ™ 
and especially a s x c ^ lar in NUHM are enhanced a lot compared to mSUGRA scenario. 
An important point we notice from the figure is that a' X p lar is usually (especially in 
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the region of large cross section) larger than a^ n /(2m x /GeV) in the NUHM case, 
while the opposite is true for the mSUGRA case. This fact implies the capture rate 
for the sun in the NUHM is largely determined by the spin-independent scattering 
cross section rather than spin-dependent one, unlike the mSUGRA scenario. This 
is because the ratio of the contribution from spin-independent and spin-dependent 
cross section to the capture rate for the sun is approximately proportional to the 
ratio of a scalar and a%j n /(2m x /GeV). 




(a) (b) 

Figure 3: the muon fluxes from (a) the sun and (b) the earth vs. B{B S — > in 
Non-universal Higgs mass scenario. 

As we have shown in the previous paper 0, the current experimental limit of 
B(B S — > puts a strong constraint on the allowed range of the spin-independent 

cross section. Since the muon flux from the sun and the earth strongly depends on 
the spin-independent cross section, we naturally expect that the current limit of 
B(B S — > play an important part in restricting the muon flux. This point 

can be observed clearly in Fig. [| where we show explicitly the correlation between 
B(B S — > and the muon flux from the sun (a) and the earth (b) in NUHM 

for tan/5 = 35 and 50. Note that the B s — > fi + fi~ is stronger for larger tan/3, 
and the resulting muon flux becomes smaller for the larger tan/3 case, like the spin 
independent DM scattering cross section @. 

The enhancements of the neutralino DM scattering cross sections (both spin- 
dependent and spin- independent) lead to the substantial change of the muon flux 
both from the sun and the earth compared with the mSUGRA case. Fig.s f| (a) and 
(b) show the muon fluxes from the sun and the earth, respectively, in non-universal 
Higgs mass scenario with 5u d — — 1, Sh u = +1 for tan (3 = 35 case. Now the maximal 
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(c) (d) 

Figure 4: the muon flux from the sun and the earth vs. m x in Non-universal Higgs mass 
scenarios with 5n d = — 1, &h u = +1 and tan/3 = 35 and tan/3 = 50. The red points (the 
open circles) are exclued by the current upper limit of B{B S — > 



values of the muon fluxes from the sun and the earth are ~ 10 3 (10) km~ 2 yr~ l , 
which is two (eight) orders of magnitude lager than the one for the mSUGRA case 
with tan/3 = 35. In Fig. ^ (c) and (d), we show the muon flux from the sun and the 
earth, respectively, in non-universal Higgs mass scenario with 5u d = — 1,<5h u = +1 
for tan/3 = 50. The red points (the open circles) are exclued by the current upper 
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limit of B(B S — > ). 
3.3 D-brane model 

Next, we consider a specific D brane model where the SM gauge groups and 3 
generations live on different Dp branes [29|. In this model, scalar fermion masses 
are not completely universal and gaugino mass unification can be relaxed. Also the 
string scale is around fO 12 GeV (the intermediate scale) rather than GUT scale. 

Since there are now three moduli (Tj) and one dilaton superfields in this case, 
we use the following parametrization that is appropriate for several T moduli: 



F s 
F i 



V3 (S + S*) m 3 / 2 sin6>, 
(Tj + T*) m 3 / 2 cos 9 9j 



(3.3) 



where 9 and Oj (i = 1,2,3) with X]J®«| 2 = 1 parametrize the directions of the 
goldstinos in the S, Tj field space. Then, the gaugino masses are given by 

M 3 = v3m 3 / 2 sin 9, 
M 2 = v3m 3 / 2 0i cos 0, 

2sin6> 



Y 



V3m 3/ 2a Y {M I 



26 3 cos# Gicos^ 
ai(M/) + a 2 (M 7 ) + 3a 3 (M 7 ) 



where 



+ 



+ 



(3.4) 



(3.5) 



a y (M 7 ) ai(JW» a 2 (Mj) 3« 3 (M 7 )' 

The string scale Mj is determined to be Mj = 10 12 (5 x 10 14 ) GeV from the U(l)i 
gauge coupling ax(Mi) = 0.1(1) P^[ . Note that the gaugino masses are non universal 
in a natural way in this scenario, unlike other scenarios studied in the previous 
subsections. 

The soft masses for the sfermions and Higgs fields are given by 
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(3.6) 
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Note that the scalar mass universality in the sfermion masses and Higgs masses is 
achieved when 

sin 2 9 = - and 2 = - for % = 1, 2, 3. (3.7) 

And in this case the gaugino masses becomes also universal, when we take only 
positive numbers for the solutions. For other choices of goldstino angles, the scalar 
and the gaugino masses become nonuniversal, and there could be larger or smaller 
flavor violation in the low energy processes as well as enhanced SUSY contributions 
to the a^ USY . 

The trilinear couplings are given by 

A u = — m 3/2 [(0 2 - ©l - ©3) cos6» - sin 9] , 

A d = — m 3/2 [(0 3 - ©i - © 2 ) cos 8 - sin9] , 

A e = 0. (3.8) 

Therefore the D brane model considered in this work is specified by following six 
parameters : 

m 3/2 , tan/3, 9, ©i=i, 2 , sign(//). 

Due to the departure from the universlity of scalar masses and the proportionality 
of trilinear couplings, the flavor violation could be different from the mSUGRA case. 
For example, it is possible to have smaller b — > s transition due to the smaller ti — tn 
mixing and larger stop masses in this D— brane scenarios, so that the B(B S — > 
/i + /i~) constraint can be relaxed. This can be seen in Fig. |5| (a) and (b), where 
the large flux signals are excluded by Super-K and AMANDA II, but not by the 
B(B S — > /i + /U _ ) constraint. In this limited parameter space, one can have a large 
DM scattering cross section and the upward-going muon flux without conflict with 
the B — > fi + fi~ branching ratio. Also there is no strong correlations among these 
observables. Therefore the indirect search for the DM annihilation is complementary 
to the B s — > branching ratio in the D— brane scenarios. 



4. Conclusions 

In this work, we considered the indirect detection of the DM through the upward- 
going muon flux from the DM annihilation at the core of the sun or the earth, along 
with the upper bound on the branching ratio for the B s fi + n~ decay, in some 
general supergravity scenarios where the upward-going muon flux could be enhanced 
very much compared to the mSUGRA case. In general supergravity scenario with 
non-universal Higgs model, we found the following: 

• Both <J^p n and a s ^ laT can be enhanced a lot compared to the mSUGRA sce- 
nario, but the enhancement in the spin-independent part is much greater. 
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Figure 5: the muon flux from the sun and the earth vs. m x in a D-brane model with 
tan/3 = 50. The red points (the open circles) are exclued by the current upper limit of 
B(B S -» fi+fi-). 



• Therefore, contrary to the usual claim, the upward-going muon flux from the 
sun can be dominated by the spin-independent part a s ^ lar in the NUHM, rahter 
than by the spin-dependent part cr^P™, as in the mSUGRA scenario [ Fig. 2 (a) 
and (b) ]. 

• The current upper bound B(B S — > fi + fi~) < 4.1 x 10" 7 excludes a large pa- 
rameter space where the muon fluxes could be enhanced otherwise, and the 
constraint is stronger for larger tan/3 [ Fig. 3 (a) and (b) ]. 

• The upper bound on B(B S — > becomes much stronger than the upper 
limits on the muon flux from Super-K and AMANDA II [ Fig. 4 (a)-(d) ]. 

In the D— brane models with nonuniversal scalar fermion masses, the correlations 
between the muon flux and B(B S — > / u + ^) becomes lost, and the upper bound on 
B(B S — ► is complementary to the upper bounds on the muon fluxes from 

Super-K and AMANDA II. Our study shows that the muon flux originated from the 
DM annihilation in the sun could be in the range of a few xlO 3 /km 2 - yr. 

Our study indicates that it is most important to include the B s — > fi + pT branch- 
ing ratio constraint when we study the direct and the indirect detections of the neu- 
tralino DM in general supergravity scenarios. The upper limit on the B s — > 
branching ratio excludes significant part of parameter space where the DM scattering 
cross section and the upward-going muon flux could be enhanced above/around the 



current experiments. Unless the chargino-stop contribution to B s — > is very 

small or there is fortuitous cancellation between the chargino-stop and the gluino- 
sbottom loop contributions, the spin-independent DM scattering cross section and 
the indirect detection rate through the upward-going muon flux are strongly con- 
trained by the B s — > branching ratio. Since both the direct and the indirect 
detection rates are well below the current experiments in most supergravity model 
parameter space when the B s — > ji + ji~ branching ratio constraint is imposed, it 
would be a great challenge for experimemtalists to reach such sensitivity to have 
positive signals of the DM search. 
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